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ABSTRACT 

In spite of recent detections of magnetic fields in a number of 8 Cephei and slowly 
pulsating B (SPB) stars, their impact on stellar rotation, pulsations, and element dif- 
fusion is not sufficiently studied yet. The reason for this is the lack of knowledge of 
rotation periods, the magnetic field strength distribution and temporal variability, and 
the field geometry. New longitudinal field measurements of four B Cephei and candi- 
date B Cephei stars, and two SPB stars were acquired with FORS 2 at the VLT. These 
measurements allowed us to carry out a search for rotation periods and to constrain the 
magnetic field geometry for four stars in our sample. 

Subject headings: stars: early-type — stars: magnetic field — stars: oscillations 
— stars: variables: general — stars: fundamental parameters — stars: individual 
CMa, 15CMa, a Pyx, eLup, 33Eri, HYVel) 



'Based on observations obtained at the European Southern Observatory (ESO programme 084.D-0230(A)). 
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1. Introduction 



For several years, a magnetic field survey of main-sequence pulsating B-type stars, namely 
the slowly pulsating B (SPB) stars and /3 Cephei stars, has been undertaken by our team with 
FORS 1 in spectropolarimetric mode at the VLT, allowing us to detect in four ft Cephei stars and 
in 16 slowly pulsati ng B stars, for the first tim e, longitudinal magnetic fields of the order of a 
few hundred Gauss (IHubrig et al.l 120061 12009|) . J3 Cephei variables have spectral types B0-B2 
and pulsate in low-order pressure and gravity modes with periods between 2 and 6 hours. Slowly 
pulsating B (SPB) stars are mid-B type (B3-B9) objects pulsating in high-order gravity modes with 
periods in the range of 0.5-3 d ays. Pulsating stars are c urrently considered as promising targets 
for asteroseismic analysis (e.g. IShibahashi & Aertsll2000|) . which requires as input the observed 
parameters of the magnetic field topology. Earl y magnetic field searches of B 'Ceph ei stars were 
mostly unsuccessful due to low precision (see lBabcocklll958l : iRudy & Kemplll978l) . Before we 
started our systematic search for magnetic fields in pulsating B-type stars, a weak magnetic field 



was d etected in two 8 Cephei stars, in the protot ype of the class, /?Cep itself, by iHenrichs et al. 



J2000l) and in V2052 Oph by iNeiner et all d2Q03ah The firs t detection of a weak magnetic field 
in the SPB star £Cas was reported by INeiner et al.1 (|2003bl) . The detected magnetic objects for 
which we gathered several magnetic field measurements showed a field that varies in time, but no 
periodicity could be derived yet due to the limited amount of VLT observing time. Among these 
targets with a detected magnetic field, we selected two slowly pulsating B stars, two f3 Cephei 
stars, and two candidate ft Cephei stars with suitable coordinates, for successive VLT multi-epoch 
magnetic measurements. The list of the selected targets is presented in TableQ] In the four columns 
we list the HD number, another identifier, the spectral type retrieved from the SIMBAD database, 
the pulsating type, and membership in a spectro scopic binary system. An asterisk in front of the 
HD number denotes candidate fi Cephei stars (cf. Stankov & Handlerll2005 ). This most recent study 
aimed at the determination of magnetic field properties for these stars, such as field strength, field 
geometry, and time variability. Here we present the results of 62 new magnetic field measurements 



Table 1: The observed fi Cephei and SPB stars. 



HD 


Other 


Spectral 


Comments 




Identifier 


Type 




24587 


33Eri 


B5V 


SPB, SB1 


46328 


g CMa 


B1III 


ySCep 


50707 


15CMa 


Bllb 


/3Cep 


* 74575 


or Pyx 


B1.5III 


/3Cep 


74560 


HYVel 


B3IV 


SPB, SB 1 


* 136504 


eLup 


B2IV-V 


ySCep, SB 
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of the six selected stars and discuss the obtained results on their rotation periods and magnetic field 
geometry. 



2. Magnetic field measurements and period determination 



Multi-epoch time series of polarimetric spectra of the pulsating stars were obtained with 
FORS2| on Antu (UT1) from 2009 September to 2010 March in service mode. Using a slit width 
of 0/4, the achieved spectral resolving power of FORS 2 obtained with the GRISM 600B was 
about 2000. A detailed description of the assessment of the longitudinal magnet ic field measure- 
ments using FORS 2 is presented in our previous papers (e.g., Hu brig et al. and references 
therein). The mean longitudinal magnetic field, (B z ), was derived using 



V 

7 



geseA 2 1 d/ 
4nm„c 2 I dA 



(1) 



where V is the Stokes parameter which measures the circular polarisation, / is the intensity in the 
unpolarised spectrum, g e g is the effective Lande factor, e is the electron charge, A is the wavelength, 
m e the electron mass, c the speed of light, dl/dA is the derivative of Stokes /, and (B z ) is the mean 
longitudinal magnetic field. The measurements of the longitudinal magnetic field were carried out 
in two ways, using the whole spectrum «5 z ) all ) and using only the hydrogen lines «5 z ) hyd ). 

Two additional polarimetric spectra of ^ CMa were obtained with the SOFIN spectrograph 
installed at the 2.56 m Nordic Optical Te lescope on La Palma , one on 2008 September 13, and an- 
other one on 2010 January 01. SOFIN (ITuominen et al.lll999|) is a high-resolution echelle spectro- 
graph mounted at the Cassegrain focus of the NOT. The star was observed with the low-resolution 
camera with R = A/AA ~ 30 000. We used the 2K Loral CCD detector to register 40 echelle orders 
partially covering the range from 3500 to 10 000 A with a length of the spectral orders of about 
140 A. Two such exposures with quarter- wave plate angles separated by 90° are necess ary to deriv e 
circularly polarised spectra. The spectra were reduced with the 4A software package (|Ilyinl l2000). 



A frequency analysis was performed on the longitudinal magnetic field measurements (B 7 ) a » 
(which generally show smaller sigmas) available from our previous (|Hubrig et al.1 12006L 120091) 
and the current studies using a non-linear least-s quares fit of the multiple harmonics utilizing the 
Levenberg-Marquardt method (|Press et al.lll992|) with an optional possibility of pre-whitening the 
trial harmonics. To detect the most probable period, we calculated the frequency spectrum for the 
same harmonic with a number of trial frequencies by solving the linear least-squares problem. At 



2 The spectropolarimetric capabilities of FORS 1 were moved to FORS 2 in 2009. 
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Table 2: Magnetic field measurements of /JCephei stars with FORS 1/2 and SOFIN (marked with 
an asterisk). All quoted errors are l<x uncertainties. 



MJD 


Rotation <Bz>„ll 






Phase [G] 


[G] 


f J CMa 



53475.046 


0.702 




282+42 


280+44 


53506.971 


0.351 




278+43 


330+45 


54061.325 


0.715 




287+42 


360+45 


54107.266 


0.795 




312+43 


319+46 


54114.028 


0.898 




309+35 


347+38 


54114.182 


0.969 




364+35 


382+47 


54116.108 


0.853 




307+45 


276+58 


54155.086 


0.738 




308+47 


349+35 


54343.371 


0.132 




345+11 


379+15 


54345.338 


0.034 




366+11 


400+12 


54345.414 


0.069 




340+11 


378+18 


54548.982 


0.476 




277+55 


297+87 


54549.995 


0.941 




380+37 


332+55 


•54722.274 


0.991 




386+39 




55107.342 


0.678 




229+30 


302+44 


55109.325 


0.589 




206+31 


233+52 


55113.224 


0.378 




203+44 


320+65 


55135.200 


0.461 




213+39 


240+59 


55150.342 


0.409 




176+51 


322+76 


55153.340 


0.784 




295+61 


470+94 


55159.329 


0.532 




207+29 


254+41 


55163.085 


0.256 




282+35 


389+56 


55164.092 


0.718 




272+45 


416+88 


55165.106 


0.183 




301+39 


431+64 


55168.091 


0.553 




232+44 


174+59 


•55201.279 


0.781 




297+26 




aPyx 


54082.341 


0.816 




142+48 


219+60 


54109.150 


0.200 




132+50 


184+60 


55107.378 


0.362 




5+31 


-29+48 


55118.347 


0.792 




120+28 


89+34 


55120.351 


0.418 




-14+26 


-35+44 


55168.167 


0.371 




24+39 


-26+68 


55171.162 


0.308 




40+45 


42+69 


eLup 


54344.998 






-156+34 


-128+36 


55225.268 






-130+94 


-191+129 


55226.272 






-33+39 


-40+52 


55227.285 






-104+44 


-185+69 


55228.324 






-7+42 


7+51 


55258.206 






-147+46 


-162+54 


55259.260 






-105+41 


-158+66 


15 CMa 


54107.318 


0.085 




163+52 


157+58 


54345.372 


0.917 




149+19 


123+27 


55107.326 


0.192 




86+26 


113+55 


55109.371 


0.354 




-41+36 


-130+52 


55112.380 


0.592 




-92+48 


-88+74 


55150.330 


0.594 




-75+54 


-49+78 


55159.295 


0.304 




35+34 


13+51 


55163.101 


0.605 




-25+39 


-45+58 


55164.104 


0.684 




31+49 


-24+61 


55165.118 


0.764 




76+41 


19+60 


55168.104 


0.000 




128+42 


126+65 


55170.090 


0.158 




138+52 


90+120 


55171.133 


0.240 




21+42 


-12+67 


55173.143 


0.399 




-58+88 


-64+122 


55177.323 


0.730 






0+62 



-5- 



Table 3: Magnetic field measurements of SPB stars with FORS 1/2. All quoted errors are la- 
uncertainties. 



MJD 


Rotation 


(#z)all 


(^z)hyd 




Phase 


[G] 


[G] 


33Eri 


52971.071 


0.599 


-122±64 


-120±68 


53574.415 


0.157 


-14+33 


-16+36 


53630.250 


0.796 


-34+27 


-32+29 


54086.175 


0.135 


-30+54 


-45±61 


54343.301 


0.098 


67±60 


83±65 


55107.154 


0.105 


116+39 


166+49 


55108.197 


0.929 


-25+43 


-43±56 


55109.350 


0.822 


5±35 


10±49 


55110.181 


0.471 


-98+48 


-138+70 


55111.222 


0.284 


-39+45 


-52+69 


55112.347 


0.164 


0±43 


-20±72 


55113.188 


0.821 


25+44 


62+73 


55120.099 


0.223 


-63+55 


-177+79 


55135.184 


0.013 


90+52 


112+64 


55149.176 


0.948 


117+37 


119+44 


55150.302 


0.829 


12±46 


-7±72 


55161.089 


0.259 


-126+46 


-168+55 


55163.031 


0.777 


-41+44 


-18+59 


HYVel 


53002.141 




-180±57 


-199+61 


53143.986 




-48+60 


-53±66 


54108.348 




-198+55 


-191+58 


55107.361 




34+59 


-61±79 


55112.367 




-25±51 


-58±69 


55118.320 




133±62 


162+83 


55120.324 




63±33 


59+48 


55168.150 




-99±38 


-147+64 


55169.349 




47+41 


78+73 


55171.150 




40+59 


-11+88 


55173.269 




17+51 


-19+68 


55177.341 




117+35 


160+49 


55181.141 




42+36 


35+54 


55182.280 




67±38 


144±53 


55189.176 




19±44 


30±68 
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each tria l frequency we performed a statistical test of the null hypothesis for the absence of peri- 
odicity (ISebenll977|) . i.e. testing that all harmonic amplitudes are at zero. The resulting F-statistics 
can be thought of as the total sum including co vari ances of the ratio of harmonic amplitudes to their 
standard deviations, i.e. as a signal-to-noise ratio (Ilyin 2000). T he F-statistics al lows to derive the 
false alarm probability of the trial period based on the F-test (|Press et al.lll992|) . Periodicity was 
found for four out of the studied six stars. The derived ephemeris for the detected periods are 



CMa : <V&/) raax = MJD55 140.73332 + 0.03794 + 2.17937 ± 0.00012£ 

1 5 CMa : < V&/) max = MJD55 1 68 .099 11 + 0.1 6667 + 12.64115 ± 0.00822£ 

a Pyx : <V&/) max = MJD55 144.59481 + 0.04105 + 3.19779 ± 0.000 19E 

33 Eri : <V&/) max = MJD55123.65285 + 0.03243 + 1.27947 ± 0.00005£ 

The logbook of the new FORS 2 and the old, revisited, FORS 1 spectropolarimetric obser- 
vations is presented in Tables |2] and [3] In the first column we indicate the MJD value at mid 
exposure. The phases of the measurements of the magnetic field, if available, are listed in Col- 
umn 2. In Columns 3 and 4 we present the longitudinal magnetic field (5 z ) all measured using the 
whole spectrum and the longitudinal magnetic field (B z ) hyd using only the hydrogen lines. Phase 
diagrams of the data folded with the determined periods are presented in Fig. \\\ The quality of our 
fits is described by a reduced ^ 2 -value which appears in the four panels of Fig. Q] 



3. Characterisation of the magnetic field geometry for the stars with determined periods 



The determination of the fundamental p arameters and the description of the pulsating prop- 
erties of all studied stars were presented by iHubrig et al.l (|2009l Tables la, lb, and 4). The most 



Table 4: Magnetic field models for the stars with detected periods. 



Object 


f 1 CMa 


15 CMa 


a Pyx 


33 Eri 


(B z ) 


[G] 


281.6+4.3 


45.3+7.1 


90.3+5.1 


-44.1+15.3 




[G] 


80.3+5.8 


116.3+12.1 


119.7+12.8 


112.0+25.5 


v sin i 


[km s^ 1 ] 


9+2 


34+4 


11+2 


25+4 


R 


[Re] 


7.1+0.9 


10.0+1.5 


6.3+1.0 


2.5+0.3 


V eq 


[kms" 1 ] 


165+21 


40+6 


100+16 


99+12 


i 


[°] 


3.1+0.8 


58.1+17.6 


6.3+1.5 


14.6+3.0 




[°] 


79.1+2.8 


57.9+18.3 


85.2+1.3 


84.1+2.7 


B d 


[G] 


5300+1100 


570+50 


3800+500 


1500+350 
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Fig. 1. — Phase diagrams with the best sinusoidal fit for the longitudinal magnetic field measure- 
ments. The residuals (Observed - Calculated) are shown in the lower panels. The deviations are 
mostly of the same order as the error bars, and no systematic trends are obvious, which justifies a 
single sinusoid as a fit function. 
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simple modeling of the magnetic field geometry is based on the assumption that the studied stars 
are oblique dipole rotators, i.e. their magnetic field can be approximated by a dipole with the 
magnetic axis inclined to the rotation axis. 

The magnetic dipole axis tilt 8 is constrained by 



r = 



(B z ) mm cos 8 cos i - sin 6 sin i 
(5 z ) max ~ cos/Jcos/ + sin/3 sin?' 



(2) 



so that the obliquity angle 6 is given by 



6 = arctan 



1 + r 



cot i 



(3) 



In Table @] we show for each star with detected periodicity in Rows 2 and 3 the mean value 
(B z ) and the amplitude o f the fi eld variation A<s z >. In Row 4 we present v sin i values published 



recently by Lefever et all (|2010|) for % l CMa and 15C Ma and from 



and 33 Eri. The radius values in Row 5 were taken from 



Hubrig et 



(R= 10.0+ 1.5 Rq) w as derived in the sa me way as in iHubrig et all t00$) by adopting the values of 



al. 



Hub rig et al.l (120091) for a Pyx 



2009). The radius of 15 CMa 



r eff and logg from iLefever et al.1 (|2010|) . In the last four rows we list the v eq and the parameters 
of the ma gnetic field dip ole models. The polar field strength in the last row was calc ulated 
following IPrestonl (U.969) using limb darkening parameters from lDiaz-Cordoves et al.l (| 19951) . 



4. Discussion 

Using FORS 1/2 and SOFIN longitudinal magnetic field measurements collected in our re- 
cent studies, we were able to determine rotation periods and constrain the field geometry of two 
8 Cephei stars, one candidate 8 Cephei star, and one SPB star. The dipole model provides a satis- 
factory fit to the data and among the very few presently known magnetic 8 Cephei s tars, $ l CMa 



and a Pyx possess the largest magnetic fields, with a dipole strength of several kG. iBriquet et al 



d2007l) discussed the position of SPB and chemically peculiar Bp stars in the H-R diagram, indi 



eating that the group of Bp stars is significantly younger than the group of SPB stars. A similar 
conclusion was deduced by IHubrig et al. ( 2007 ). who studied the evolution of magnetic fields in 



Ap and Bp stars with definitely determined magnetic field geometries across the main sequence. 
The vast majority of Bp stars exhibits a smooth, single-wave variation of the longitudinal magnetic 
field during the stellar rotation cycle. These observations are considered as evidence f or a dominant 



dipo lar contribution to the magnetic field topology. It is of interest that the study of Hu brig et al 



(|2007|) indicates the prevalence of larger obliquities 8, namely 8>60°, in more massive stars. The 
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magnetic field models for the three /3 Cephei stars and the one SPB star presented in this work 
confirm this trend. 

The insufficient knowledge of the strength, geometry, and time variability of magnetic fields 
in hot pulsating stars prevented until now important theoretical studies on the impact of mag- 
netic fields on stellar rotation, pulsations, and element diffusion. Althoug h it is expected that the 
magnetic field can distort the frequency patterns (e.g. lHasan et al.1 120051) . such a perturbation is 
not yet d etected in hot pulsating stars. Splitting of non-radial pulsation modes was observed for 
15CMa (|shobbrook et al. 2006 ). but the identification of these modes is still pending. The mag- 
netic 6 Cephei star samp le indicates that they all share common p roperties: they are N-rich targets 
(e.g., iMorel et al]|2008|) and, as discussed by Hubrig et al. I J2009I) . their pulsations are dominated 
by a non-linear dominant radial mode (see also Saesenl 20061 for £ l CMa). The presence of a 
magnetic field might consequently play an important role to explain such a distinct behaviour of 
magnetic f3 Cephei stars. More precisely, chemical abundance anomalies are commonly believed 
to be due to radiatively-driven microscopic diffusion in stars rotating sufficiently slowly to allow 
such a process to be effective. However, we need an additional clue to account for the fact that both 
normal and nitrogen-enriched slowly rotating stars are observed. The presence of a magnetic field 
is a very plausible e xplanation, as i t can add to the stability of the atmosphere, allowing diffusion 
processes to occur (|Michaudlll970|) . On the other hand, among the studied stars, apart from the 



star 15 CMa with rather low v eq = 40 + 6 km s 1 , the other three magnetic pulsating stars rotate 
much faster up to v eq = 165 ± 21 km s _1 for £ ! CMa with the strongest magnetic field, indicating 
that these stars are not truly slowly rotating stars, but seen close to pole-on. Obviously, the topic 
of mixing signatures is not understood theoretically yet and more computational work as well as 
future additional observational validation of our results are needed to understand the link between 
the presence of a magnetic field, rotation, pulsating characteristics, and abundance peculiarities. 



T.M. acknowledges financial support from Belspo for contract PRODEX GAIA-DPAC. 
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